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Abstract 

The main component of the s process is produced by low mass stars (1.5 < M/Mq < 
3), when they chmb for the second time the red giant branch and experience a series 
of He shell flashes called thermal pulses. During the relatively long period (10^ yr) 
that elapses between two subsequent thermal pulses, a slow neutron flux is provided 
by the ^^C(a,n)^^0 reaction taking place within a thin ^^C pocket located in the 
He-rich and C-rich mantel of these stars. A second, marginal, neutron burst occurs 
during the thermal pulse and it is powered by the ^^Ne(a,n)^^Mg reaction. We 
review the present status of the nucleosynthesis models of low mass AGB stars. The 
advance in the knowledge of the complex coupling between convective mixing and 
nuclear process, which allows the surface enrichment of C and s-process elements, is 
presented, together with the hypotheses concerning the physical mechanism driving 
the formation of the ^^C pocket. In order to illustrate the capabilities and the 
limits of the theory, an updated computation of a 2 M© stellar structure with solar 
chemical composition is reported. This model has been obtained by including a full 
nuclear network (from H up to Bi, at the termination point of the s-process path) 
into the stellar evolution code. The predicted modiflcation of the surface composition 
occurring during the AGB evolution is shown. The new challenge of AGB modeling, 
namely the study of C-rich and s-rich very metal-poor stars, is discussed. 
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1 Introduction 



It was in 1868 wlien tlie Jesuit astronomer Father Angelo Secchi first rec- 
ognized in th e pec uliar spectrum of some red giants the signature of carbon 



enhancement This new class includes stars belonging to the Disk of the 

Milky Way, called by Secchi "red carbon stars" and nowaday named C(N- 
Type) stars, that are evolved low mass red giants with photospheric C/0 > 
1. In the Hertzsprung- Russell diagram, they are located near the tip of the 
asymptotic giant branch (AGB) and represent the end point of the evolution- 
ary sequence of a star that starts the AGB as an M giant (typically C/0 ~ 
0.5) and progressively modifies its surface composition, passing through the 
MS, S and C(N) stages. The internal structure of these giant stars is made 
of three regions: a compact C-0 core, a thin He-rich and C-rich mantel (He 
intershell) and an expanded H-rich envelope. As usual for red giant stars, the 
envelope is largely unstable against convection. Carbon is synthesized by the 
3a reactions that burn at the base of an He-rich layer surrounding the core. 
The first question for the theoreticians was the search for a process capable 
tomove the C from the deep interior of the star to the surface (jo^, |87| . 
(ll3l |). Modern studies have clarified that the dredge up of C is due to the 
combination of two distinct convective episodes. The first is responsible for an 
efficient mixing of the whole He-rich layer and the second partially overlaps 
the zone previously mixed by the first convective episode and extends to the 
stellar surface. As a consequence of these two convective episodes, other prod- 
ucts of the nucleosynthesis occurring within the He-rich mantel, besides ^^C, 
should appear at the surface. As a matter of fact, MS, S, C(N) and post- AGB 
stars are enriched in s-process elements, like Sr, Y, Zr, and Ba, La, Ce, Nd 
corresponding to the light and heavy s-abundance peaks, as early encoded by 
[l6| . Detection of unstable ^^Tc, whose half life is "only" 2.1 x 10^ yr, demon- 
strates that this enhancement cannot be due to an anomalous pollution of 
the pristine material from which these stars were born, but that the s process 
must be at work in their interiors In addition, it was realized that only 
relative ly fa int AGB stars become C-rich and that the brightest AGBs are 
N-rich |l29j . This evidence supports the hypothesis that C(N)-stars have low 
mass progenitors. In massive AGBs, indeed, the CN-cycle taking place at the 
base of the external convective layer converts most of the C dredged up into 
N. 

Recent studies have demonstrated that the AGB stars with M < 3 Mq are 
the major contributors to the galactic production of the isotopes belonging 
to the main component of the s p roce ss, involvin g iso topes with atomic mass 



beyond A ~ 80 (see for a review [128[). Clayton j31| first pointed out that a 



multiple neutron exposure, rat her t han a single exposure, is needed to repro- 



duce the main component and |l2lj suggested that this condition could occur 



in the He-rich layer of an AGB star. A few years before, indeed, 110l | discov- 
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ered that during the late AGB the He burning shell is recurrently switched 
on and off. At each He reignition, a thermonuclear runaway (thermal pulse) 
occurs and the rapid release of nuclear energy induces the formation of an ex- 
tended convective region, where the products of the nucleosynthesis are fully 



mixed. Following the original idea of 



121 



the partial overlap of the convec- 



tive regions recurrently generated by the thermal pulses ensures the required 
multiple neutron exposure (see Figure 1). 
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Fig. 1. This sketch illustrates the evolution of the positions of the inner border of 
the convective envelope, the H-burning shell and the He-burning shell, during the 
thermally pulsing AGB phase. The convective regions generated by two subsequent 
thermal pulses are also shown. Note that the temporal developments of thermal 
pulses and of following TDU episodes are off scale with respect to the interpulse 
period (see Figure 3 and 8). 



Two promising neutron sources for the build up of the heavy s-process isotopes 
were early recognized (t21l], 12^]): the ^^Ne(a,n)^^Mg and the ^^C(a,n)^^0 re- 
actions. On the basis of detailed AGB stellar models, it has been understood 
that the major source of neutrons in low mass AGB stars is provided by the 
^^C(a,n)^^0 reaction. The s-process nucleosynthesis mostly occurs during the 
relatively long interpulse period (namely the time elapsed between two subse- 
quent thermal pulses), in a thin radiative layer at the top of the He intershell 
(~ 10~^ ^o); when the temperature ranges between 80 and 100 x 10^ K [11^. 
A second neutron burst giving rise to a small neutron exposure ^ , but with a 
high peak neutron density, is released by the marginal activation of the ^^Ne 

^ We define the neutron exposure r = / rinVthdt as the time integrated neutron 
flux, where n„ is the neutron density and Vth the thermal velocity. 
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neutron source in the convective thermal pulse, modifying the final s-process 
composition at branchings along the s path depending on the neutron density 
or on temperature. 

In this paper we discuss the present status of the extant studies of nucleosyn- 
thesis and evolution in low mass AGB stars. The AGB evolution is reviewed 
in sections 2 and 3. In section 4 the general problem of the third dredge up 
efficiency in low mass AGBs of different mass and metallicity is discussed. An 
attempt to properly evaluate the mass loss rate, a quantity that significantly 
affects the AGB evolution and nucleosynthesis, is presented in section 5. In 
section 6 the problem of the source of neutrons is afforded in some detail, 
with particular emphasis to the physical process driving the formation of a 
^■^C pocket in the He-rich and ^^C-rich intershell. An AGB stellar model of 
initial mass M = 2 Mq and solar metallicity, computed by coupling a stel- 
lar evolution code with a full nuclear network (from H to Bi), is described in 
section 9. Finally, among the many aspects related to the s-process nucleosyn- 
thesis occurring in AGB stars, we focus, in section 10, on one of the presently 
most debated problems in the field, i.e. the expected s-process nucleosynthesis 
for the most metal-poor stars in the Milky Way. 



2 Approaching the Thermally Pulsing AGB 



In stars less massive than ~ 10 Mq, central He burning leaves a compact 
C-0 core. After He exhaustion, the central density rapidly increases (above 
10^ and up to 10^ g/cm~^), electrons become highly degenerate and a huge 
energy loss by plasma neutrinos takes place. The neutrino energy depletion 
is only partially supplied by the release of gravitational energy from the core 
contraction, and the thermal content of the core is used to balance the deficit. 
The resulting cooling starts from the center, where the density is higher, and 
the maximum temperature moves progressively outward. It exists a maximum 
mass, called M„p, for which the whole core cools down, a fact that prevents 
carbon ignition. The precise value of M^p depends on the chemical composi- 
tion: it is ~7 Mq for population I stars (of nearly solar composition) and for 
very metal-poor stars (population HI), while it is smaller (~6 Mq) for popu- 
lation II (halo) stars. Stars with mass slightly large r than Mup suffer a violent 
carbon ignition in degenerate condition (|40| Jl03|). Stars with smaller mass 
enter the thermally pulsing AGB phase (see [38[ for an updated presentation 
of the pre- AGB evolution). 

The evolutionary track in the HR diagram calculated for a 2 Mq star of solar 
metallicity is reported in Figure 2. The computation has been started from a 
homogeneous structure with a mass fraction of helium Y = 0.275 and a metal- 
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licity Z = 0.015^. The mass fractions of all the elements beyond hydrogen 
and helium relative to Z have been derived from 0], except for C and N 
and O, Ne and Ar 0]. The initial model is fully convective and corresponds to 
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Fig. 2. Theoretical evolutionary track in the Hertzsprung- Russell diagram of a star 
of initial mass M = 2 Mq and solar metallicity (see section 9 for details of the 
model). 

the pre-main sequence contraction phase, with a central temperature of about 
10^ K. The model has been evolved from the pre-main sequence, through core 
H burning, the red giant branch (RGB), the off-center He flash and the He 
burning, up to the AGB. 

After the star leaves the main sequence and first becomes a red giant, the con- 
vective envelope penetrates in the radiative region above the H shell, enriching 



^ These values of Y and Z have been derived from a standard solar model, by 
fitting the luminosity (Lq = 3.844 • 10^^ erg/s), radius {Rq = 6.951 • lO^'' cm) and 
the mass fraction of metals relative to hydrogen {{Z/X)q = 0.017) of the present 
Sun. 
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the surface composition with the ashes of proton captures that occurred in this 
deep zone during the main sequence (first dredge up, FDU); with respect to 
the pristine composition, ^He, '^He, ^^C, ^JO, and ^^N are enhanced, while 
i2q^ i5]\j g^j^^ 18Q depleted (see e.g. |l5|). Light elements, Li and Be, are 
practically extinguished. The upper mass limit for the occurrence of the first 
dredge up increases with the metallicity. At Z = 0.02, all stars entering the 
AGB have experienced the first dredge up episode, but this only occurs for 
M < 3 Mq and M < 2.5 M© when Z = 0.001 and 0.0001, respectively ([ij, 
0). 

During the first part of the asymptotic giant branch (called early- AGB) , the 
He shell burning progressively moves outward and the mass of the C-0 core 
increases. In low mass stars, the H burning maintains an entropy barrier that 
limits the internal boundary of the external convective layer. In contrast, in 
massive AGB stars, owing to the huge energy flux coming out from the He 
burning zone, the base of the H-rich envelope expands and cools, so that the 
H burning dies down. In this case, the external convection penetrates inward, 
within the H-depleted zone. This is the second dredge up (SDU), found in 
stellar models with M > 3—5 Mq, depending on the chemical composition^ . 
As a consequence of the SDU, a further increase of the surface abundances 
of helium and nitrogen is expected in massive AGBs. In addition, the SDU 
reduces the H-depleted region and prevents the formation of massive white 
dwarfs. 

Finally, when the He burning shell gets closer to the H/He discontinuity, it 
dies down and, after a rapid contraction, the H burning shell fully supplies 
the surface energy loss. 



3 Thermally pulsing AGB 



The temporary stop of the He burning shell marks the beginning of the ther- 
mally pulsing AGB phase (TP- AGB). The first thermal pulse occurs when 
the H burning shell accumulates enough He below it, so that the He-rich zone 
is compressed and heated, and He reignites. Although the degree of electron 
degeneracy of the He-rich material is weak, a thermonuclear runaway occurs, 
because the thermodynamical time scale needed to locally expand the gas is 
much longer than the nuclear burning time scale of the 3a reaction jl09l |. Ow- 
ing to the fast release of nuclear energy, the local temperature increases and 
the He burning luminosity blows up, in extreme cases to 10^ Lq. The ther- 
monuclear runaway drives the formation of a convective zone that extends 

^ the lower the metallicity, the smaller the minimum mass for the occurrence of the 
second dredge up. 
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from the region of the partial He burning to the H/He discontinuity (see Fig- 
ure 3). At the base of this convective shell an incomplete He burning takes 
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Fig. 3. Convective episode (dashed area) in the model of M = 2 Mq, Z = 0.015 
and Y = 0.275, during and after the lO*'^ thermal pulse. Note how the convective 
zone generated by the thermal pulse covers the whole He intershell. After about 200 
years, the external convection penetrates inward (TDU). 

place and the products of the 3a reaction (essentially carbon) are mixed over 
the whole intershell. At the quenching of a thermal instability, the resulting 
mass fraction of C in the top layer of the He intershell is X(^^C) ~ 0.25. When 
the expansion has progressed far enough, the temperature of the He shell de- 
creases and a quiescent He burning phase begins. The variations of the H and 
He burning luminosities for the 2 Mq model with Z = 0.015 are reported in 
Figure 4. The expansion powered by the thermal instability pushes outward 
the layers of material located above the He burning shell. The temperature 
and the density at the base of the H-rich envelope decrease and the H burn- 
ing shell dies down. In these conditions, a third dredge up (TDU) episode is 
strongly favoured. Indeed, at the interface between the envelope and the man- 
tel, owing to the large energy flux coming from below, the local temperature 
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Fig. 4. Thermally pulsing AGB evolution of the model M = 2 Mq, Z = 0.015 and 
Y = 0.275; panel a): energy production rate of the H burning shell; panel b): energy 
production rate of the He burning shell. 

gradient increases. On the other hand, due to the low density, the ratio of the 
gas pressure to the radiation pressure decreases and the adiabatic tempera- 
ture gradient approaches its minimum allowed value for a fully ionized gas 
plus radiation (i.e. {d\ogT/d\ogP)ad = 0.25). Then the Schwarzschild criterion 
for the convective instability is more easily fulfilled and the envelope may pen- 
etrate (in mass) within the He intershell. The propagation of the convective 
instability is self-sustained due to the increase of the local opacity that occurs 
because fresh H (high opacity) is brought by convection into the He-rich layers 
(low opacity). If the dredge up is deep enough to overlap the region previously 
mixed by the convective zone generated by the TP, helium, carbon and heavy 
s elements are brought to the surface. A C(N) star may event ually form as a 
consequence of a series of recurrent TDU episodes (j67|. [s^, |ll3l | ). 

Between two subsequent pulses (interpulse phase), the H burning shell sup- 
plies the energy radiated by the stellar surface and the luminosity of the star 
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basically depends on the mass of the H-exhausted core, Mh. In other words, 
a direct correlation exists between the core mass and the stellar luminosity 



9J|, [69[) • As the H burning shell advances in mass, the internal boundary 
of the convective envelope recedes, but always remains very close to the region 
where the thermonuclear reactions are switched on. 

If the stellar mass is sufficiently large (M > 5 Mq) ^ , some nuclear processes 
are indeed active within the most internal layers of the convective envelope. In 
the more massive AGBs (~ 7 Mq), the temperature at the base of the convec- 
tive envelope ma y rea ch 80 x 10^ K. This phenomenon is called hot bottom 
burning (HBB) ( |ll8l |. ji^, The main consequences of the coupling of 

nuclear burning and external convection are: i) a substantial increase of the 
stellar luminosity, which deviates from the classical core mass-luminosity re- 
lation (O, Q); ii) a modification of the surface composition. In particular. 



most of the C dredged up is converted into N, a fact that prevents the forma- 
tion of massive C(N)-stars. Moreover, the surface ^^C/^^C ratio app ro ache s 
its equilibrium value (3.5) and an enhancement of Li is possible 104 1. 123l 
through the so-called Cameron-Fowler mechanism 



4 The third dredge up 



The products of the s-process nucleosynthesis as well as those of the partial 
He burning can be actually observed only if the third dredge up takes place. 
The TDU is driven by the expansion of the envelope powered by a thermal 
pulse. It is deeper when the strength of the pulse, as measured by the maximum 
luminosity attained by the He burning, is stronger. In principle, the strength of 
the pulse depends on temperature and density of the He-rich layer at the epoch 
of the ignition. The higher is the density and the temperature at He ignition, 
the stronger the thermal pulse is. For this reason, since He is accumulated 
by the thermonuclear fusion occurring in the H burning shell at the base of 
the envelope, the H burning rate is among the most important quantities 
regulating the physical conditions of the point where He ignites. It follows 
that the third dredge up is infiuenced by the parameters affecting the H- 
burning rate such as, in particular, the metallicity, the mass of the H-exhausted 
core and the mass of the envelope H, |ii3)- As a general rule, a 

slower H burning implies a higher density of the He-rich layer and, therefore, 
a stronger thermal pulse ensues. Such an occurrence suggests a few important 
considerations: 



^ The classical core mass-luminosity relation, as reported by (ssl) relates the max- 
imum luminosity reached toward the end of the interpulse phase and the mass of 
the H-exhausted core, namely: L^ax = 5.925 • 10'' (M/^ -0.495). 
^ This lower limit is smaller at low metallicity. 
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i) along the AGB, the mass of H-depleted material that is dredged up in a 
single episode (5Mtdu) initially increases, because the core mass increases, 
reaches a maximum and then decreases, when the mass loss erodes a substan- 
tial fraction of the envelope (see Figure 9); 

ii) for given core and envelope masses, the TDU is deeper in low metallicity 
stars, because the H burning is less efficient. We found that SM^bv [Z = 
0.002) ~ 2 X SMtdv {Z = 0.02), for the same core and envelope mass. The 
relation of the TDU with the metallicity should be carefully considered when 
comparing disk AGB stars with metal-poor extragalactic AGB stars, like those 
belonging to the Magellanic Clouds or to the Dwarfs Spheroidals galaxies js^]; 

iii) there is a minimum envelope mass for which the TDU takes place. This 
minimum depends on the core mass and on the chemical composition of the 
envelope. We found that the TD U ceases when the envelope mass becomes 
smaller than 0.3 0.5 Mq 115l |. This implies that in stars of initial mass 



below a given threshold (M < 1.2 Mq), the residual envelope mass at the 
beginning of the thermally pulsing AGB is already too small and the TDU 
cannot takes place. As a matter of fact, AGB stars belonging to the Galactic 
Globular Clusters, whose initial mass are of the order of 0.8 — 0.9 Mq, do 
not show the enhancement of C and s-elements, which is the signature of the 
TDU [II2I. Similarly, at variance with their metal-rich disk counterparts, halo 
post-AGB stars do not show any significant enhancement of the s-elements 



5 Mass loss 



During the AGB, the star may become unstable against large amplitude pul- 
sations. Pulsations induce a compression of the gas and the resulting increase 
of the density of the cool atmospheric layers favours the formation of complex 
molecules and dust grains, which trap the outgoing radiation fiux driving a 
strong wind. Mass loss could also be infiuenced by the environment, as in 
close binary systems or in crowded stellar populations, like the central re- 
gion of Globular Clusters. Mass loss erodes the envelope causing important 
changes in the stellar properties. The duration of the AGB, the strength of 
the pulse, the efficiency of the third dredge up are a few examples of the quan- 
tities affected by the mass loss. The correct evaluation of the mass loss rate 
is also required to estimate the degree of chemical pollution of the interstellar 
medium ascribed to AGB stars. 

The AGB mass loss rate may be estimated from infrared color indices or from 
molecular CO rotational lines measurements. The available data indicate that 
it ranges between 10~^ and 10~^ Mq/jt. In variable AGB stars, the larger 
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the period the larger the mass loss is. No other correlations between mass 
loss and stellar parameters (luminosity, mass, composition) have been clearly 
identified (see e.g. |133l] ). The poor knowledge of the actual mass loss rate 
is one of the major uncertainties of AGB stellar models. Much effort has 
been devoted to derive a suitable prescription for this phenomenon to be 
used in calculations of AGB evolutionary models. The Reimers' form ula ^ 
was introduced to describe the mass loss in population II red giants |lOCl| 
and the rj parameter was calibrated according to the lumi nosity and color 



distribution of bright globular clusters stars {rj ~ 0.4, |101| ). Unfortunately 
an equally stringent constraint for the calibration of the mass loss rate in AGB 
stars is lacking. In principle, the mass loss rate may be adjusted in order to 
reproduce the observed luminosity functions of the AGB stellar population, 
or that of a sub-sample of AGBs, like the C(N)-stars. Another important 
constraint can be derived from the initial to final mass relation [132]. On 
the base of synthetic AGB models, 5^ suggested that a Reimers' mass loss 



formula, with the parameter rj = 5, provides a suitable reproduction of these 
observational constraints. Their method, however, only indicates the average 
mass loss rate, whereas the mass loss history remains largely unknown. Indeed 
[s?! showed how different mass loss prescriptions may equally fulfil the same 
constraints. Our tests, which make use of detailed stellar models to reproduce 
similar observational constraints, show that the mass loss rate suggested by 
[5^ may be adequate only for the more massive and/or more evolved AGB 
stars, but in low mass AGBs it is definitely too high. 

A possible alternative method to estimate the mass loss rat e is b ased on the 
observed correlation with the pulsational period (see e.g. |l22l |). Since the 
evolution of the pulsational period depends on the variations of radius, lu- 
minosity and mass, this relation provides a simple method to estimate the 
evolution of the total stellar mass from basic stellar parameters. In Figure 5 
we have collected data for the mass loss rate versus period, as measured by 
various authors for 0-rich (open circles) and C-rich (triangles and squares) 
AGB stars. Beside the evident spread of the data, three different regimes can 
be recognized: i) P < 300 days: moderate mass loss (roughly between 10~^ and 
3 X 10"'' Mq/jt); ii) 300 < P < 1000 days: exponential increase of the mass 
loss; iii) P > 1000 days: the mass loss approaches a maximum (~ 5 x 10~^ 
Mf7i/y r), w hich roughly coincides with the expected radiation-pressure-driven 
limit 122l |. The dot -das hed line in the figure represents the prescription by 



Vassiliadis & Wood 122l | . It clearly underestimates the mass loss rate of short 



period variables, whilst for P ranging between 500 and 1000 days the mass 
loss appears too high by an order of magnitude. Thus, we have worked out 
a new relation (solid line), which provides a better calibration of the mass 



3 

^ M(M0/yr) = 1.34 • 10~^ • r] ■ jff^, where L and M are in solar units and the 

cff 

effective temperature in K. 
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Fig. 5. Comparison of various mass loss rates versus period measurements (symbols) 
and prescriptions used in stellar evolution calculations (lines). 

loss-period relation, namely: 
for logP < 2.5 logM = -7.7 

for 2.5<logP<3.1 logM= -101.6 + 63.26-logP- 10.282- (logP)2 
for logP > 3.1 logM = -4.3 

This prescription has been used in our latest calculations of low mass AGB 
stellar models. In Figure 6 we report the evolution of the envelope mass (Menv) 
of an AGB model of initial mass 2 Mq, as obtained under different prescrip- 
tions for the mass loss, namely: Reimers (r/ = 5), Reimers {rj = 0.5), Vassil- 
iadis & Wood (VW) and our new cahbration of the mass loss-period relation 
(PW). The of the temporal scale corresponds to the beginning of the TP- 
AGB phase. Note how the mass loss history depends on the rate prescription. 
While the Reimers' formula provides a constant increase of the mass loss, 
the Vassiliadis and Wood rate leads to a negligible mass loss for the major 
part of the AGB evolution, with a sudden increase to the radiation-pressure- 
driven limit toward the end. Our prescription resembles a moderate Reimers 
{rj = 0.5) mass loss rate for the first 1.2 Myr and switches to the stronger 
= 5 in the late TP- AGB phase. Note how the different prescriptions affect 
the duration of the AGB and, in turn, the AGB luminosity function and the 
estimated contribution of AGB stars to the Galactic chemical evolution. 
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Fig. 6. Variation versus time of the envelope mass according to different mass loss 
prescriptions (see text), for an AGB model of initial mass M = 2 Mq. 

6 The source of neutrons and the formation of the ^'^C pocket 



As recalled in the Introduction, in thermally pulsing AGBs two major neutron 
sources may operate within the intershell, driven by the ^^Ne(a,n)^^Mg or by 
the ^^C(a,n)^^0 reaction. 

The amount of ^^N left by the H burning at the top of the intershell is prac- 
tically equal to the sum of the abundances (by number) of the C-N-0 in the 
envelope. During the early phase of the convective thermal pulse, the material 
within the intershell is fully mixed and the ^^N is totally converted into ^^Ne 
throughout the chain ^^N(a,7)^^F(/5+z/)^'^0(a,7)^^Ne. Near the peak of the 
thermonuclear runaway, if the temperature is high enough (i.e. ~3.5 x 10^ 
K), the ^^Ne(a,n)^^Mg reaction may provide a significant neutron flux. Iben 
([65[) first demonstrated that this condition is fulfilled in intermediate mass 
stars. Actually, in low mass stars (M < 3 Mq), the temperature at the base 
of the convective zone generated by the TP barely attains 3 x 10*^ K and the 
^^Ne neutron source is marginally activated (ItT 

In low mass AGB stars, an alternative neutron source is provided by the 
reaction ^'^C(Q;,n)^^0. This reaction requires a substantially lower temperature, 
namely ~90 x 10^ K, which is easily attained in the He intershell. Note that 
between two subsequent thermal pulses, the H burning shell leaves some ^'^C 
in the upper region of the He-rich intershell. The burning of this ^^C, however, 
produces a negligible neutron fiux, even if, as a consequence of the TDU, the 



nj, H). 
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CNO in the envelope may grow up to 10 times the solar abundance. This 
occurs because, in the material processed by the CNO burning, the ^^N is, in 
any case, two orders of magnitude more abundant than the ^^C. In practice, 
all the neutrons released by the ^^C left by the H burning are captured by the 
abundant ^^N. Then, an alternative source of ^^C is needed, in a zone where 
i^N is depleted. The isotope is produced via ^^C{p,-fy^N{l3+i^y^C. After 
each thermal pulse, the material within the He-rich intershell is plenty of ^^C 
and all the ^^N has been converted into ^^Ne, but it is also H depleted. So, 
the problem is how to inject a few protons into the He-rich intershell. 



Sackmann et al. lOSl ] suggested that the convective zone generated by the 



thermal pulse could extend beyond the H/He discontinuity. In this case, pro- 
tons are ingested downward, where the temperature is so high that the reaction 
chain ^^Cr2/y)^^N(/3"'"i/)^^C(a,n)^^0 would release neutrons. Against this hy- 
pothesis, 6y| argued that during the thermonuclear runaway the H burning 



shell remains active and generates an entropy barrier that prevents the pen- 
etration of the convective instability into the H-rich envelope. Extant AGB 
stellar models show that only in very metal-poor stars, owing to the l ack o f 



CNO , ingestion of protons in the convective pulse may occur (j|4^, [7l| . 
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117| ). More favorable conditions for the formation of a ^'^C pocket within 
the intershell are realized at the epoch of the TDU and during the so-called 
post-flash dip, the period that immediately follows the TDU [69|. Indeed, as 
a consequence of the TDU, a sharp discontinuity between the H-rich envelope 
and the He- and C-rich intershell forms. This condition is maintained until H 
burning is reignited. The time elapsed from the maximum penetration of the 
convective envelope and H reignition is about 10^ yr (for a star of 2 Mq). Dur- 
ing this period, different physical mechanisms may contribute to the downward 
diffusion of a few protons into the underlying radiative layer. In such a case, 
at H reignition the top layers of the He intershell heat up and a ^^C pocket 
forms. Roughly speaking, the required mixing process must be able to diffuse 
about 10^^ Mq of H within a region as large as ~ 10^'^ Mq. Note that an 
excess of protons must be avoided, because in that case the production of ^'^C 
is foll owed by the production of ^^N (via the ^^C(p, 7)^^N reaction). Straniero 
et al. m demonstrated that the -C formed in this way is fully consumed by 
the ^^C(Q;,n)^^0 reaction in radiative conditions during the interpulse phase, 
when the temperature rises up to ~ 90 x 10^ K, giving rise to a large neutron 
exposure with a maximum neutron density of about 10'' cm~^. The expected 
s-process nucleosynthesis has been first calculated by [49| . 



Among the various attempts to find the possible mechanism responsible for 
the formation of the ^'^C pocket, (6^ evaluated the timescale of the atomic 
diffusion driven by the sharp composition gradient left by the penetrating 
convective envelope at the time of the TDU and concluded that this is a 
promising possibility to diffuse enough protons from the envelope into the He- 
rich intershell. More recently [g^I, inspired by hydrodynamical simulations of 
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|45| , invoked mechanical overshoot of material from the bottom of the convec- 
tive envelope into the underlying stable region. Langer et al. Ill calculated 



AGB models taking into account stellar rotation and found a certain mixing 
of protons into the intershell at the epoch of the TDU. However, this rotation- 
ally induced mixing does not stop when the convective envelope recedes, but 
continues during the interpulse period causing the contamination of the ^'^C 
pocket with too much ^^N. Gravity waves have been investigated to explain 
the formation of the ^^C pocket by js^, and [i^l show that a magnetic field 
of 10® gauss, comparable with the one measured in white dwarfs, may induce 
a deep circulation below the convective envelope at the epoch of the TDU, 
allowing the required diffusion of protons into the intershell. 

Owing to the lack of a reliable description of the physical phenomena that 
govern the diffusion of protons into the He intershell, in the current AGB 
nucleosynthesis calculations the amount of ^^C is assumed as free parame- 
ter (see ji^). In addition, the simultaneous solution of the stellar structure 
equations and a full network including all the relevant isotopes up to the ter- 
mination point of the s-process path (Pb-Bi) requires a relevant computational 
power. For this reason, a post-process nucleosynthesis calculation, based on 
AGB stellar models computed with a restricted nuclear network, is generally 
preferred. In spite of these limitations, this approach provided a substantial 
improvement with respect to prior calculations, based on the so-called classi- 
cal analysis of the s process (see discussion in jQ]). The adequacy of detailed 
post-process calculations has been proved by comparing theoretical expecta- 
tion s wi th observed abundances found in galactic disk AGB stars ([s^, [H, 
102l |. 17 1). Additional confir mat io ns h ave been derived from presolar SiC 



grains found in meteorites ( |8q| , |10q | , |137j ) . These calculations are the brick s 
of the extant chemical evolution models incorporating heavy elements jl2(|. 
In general, observations show an important scatter of s-process abundances at 
any metallicity ^19]. This scatter reveals a spread of the effective amount of 
^^C in the pocket, perhaps due to differences in the stellar parameters, such as 
the initial mass or the mass loss rate, or to the chaotic nature of the process 
responsible of the proton injection into the He intershell. 



7 The nuclear network 



The next step in the understanding of the s-process nucleosynthesis in AGB 
stars is the development of algorithms that allow the formation of the ^^C 
pocket and a full coupling of the equations describing the physical and the 
chemical evolution of the star. In this framework, we have included a full net- 
work of about 500 isotopes (from H to the Pb-Bi-Po ending point) linked by 
more than 750 reactions, in the stellar evolution code FRANEC |22i . This net- 
work, the same already used in the post-process calculations, is continuously 
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upgraded according to the latest theoretical and experimental nuclear physics 
improvements. The adopted nuclear reactions are summarized in the following 
five subsections. 



7. 1 Charged particle capture reactions 



Reaction rates of isotopes involving charged particles are generally taken from 
the NACRE compilation j^. For those reactions not included in this database, 
we use the rates tabulated by EtI]. Two important exceptions should be re- 
called. Concerning ^^N(p,7)^^0, we use the result of the recent measurement 
obtained by the LUNA collaboration 4J|. Being the bottleneck of the CNO cy- 
cle, this reaction regulates the He deposition on the intershell of an AGB star 
and in turn affects the pulse strength and the efficiency of the TDU (see the 
discussion in section 4). The new low energy astrophysical factor (S'(0)=1.7 
MeV barn) is about 40% lower than that suggested by NACRE. With regard 
to ^^C(q;,7)^^0 we adopt the rate recently derived by [83^). 



7.2 Neutron sources 



A subthreshold resonance makes difficult the low energy extrapolation of the 
^^C(a,n)^^0 reaction rate. We adopt the rate siiggested by j4l|. We have also 



investigated the effects of different rates (|8^, [5|) and we found a neg ligible 
modification of the final surface isotopic distribution (js^], see also, [Qq). The 
^^Ne(a,n)^^Mg rate is even more complex, owing to the possible existence of 
unknown low energy resonances [t^I. We adopted the lower limit suggested 
by [l^- It (^oes not include the possible resonance at 633 keV and adopt a 
1 a lower limit (164 /leV) for the strength of the resonance at 828 keV. The 
recommended value by the NACRE compilation at T ~ 3 x 10^ K is a factor 
of two lower than our adopted choice, and very close to the recent determi- 
nation of [t^I. With such a rate, only a few neutron-rich isotopes involved in 
branchings sensitive to the peak neutron density, like *^Kr, ^"^Rb, ^^Zr, would 
be affected. 



7.3 n-capture reactions 



With n-capture reactions we mean (n, 7), (n, a) and {n.,p) processes. For the 
(n, 7] reactions, we adopt, as reference compilation, the recommended rates 
by [8] (hereafter BK2000). In that paper, experimental and theoretical data 
are critically revised and the reaction rates are listed as a function of the 
thermal energy from 5 to 100 keV. We recall that the typical thermal energy 
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at which the s-process nucleosynthesis occurs in AGB stars are 8 keV for the 
^^C(a,n)^^0 reaction and 23 keV for the ^^Ne(a,n)^^Mg reaction. For the few 
reactions not inchided in BK2000, we use the theoretical calculations of |97i]. 

Starting from this database, we have upgraded the network with the most 
recent experimental results. With regard to Si isotopes we refer to js^, CI 
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Xe isotopes to |9; 



isotopes to 



6°Ni to |33|, ^^Ni to 96], ^^Sr to 
, Cd isotopes to [135J, ^^^La to 



l81l . Kr isotopes to 
l92l |. Pm isotopes to 
and finally Pt isotopes to [7^. The rates of 



^^^Sm and Eu isotopes to |l 
{n,a) and {n,p) reactions involving heavy isotopes are taken from RT2000. 
We adopted the (n, a) and (n, p) rates on light isotopes from various authors. 
In p articular: the ^^N(n,p)^"^C is taken from js^], the ^'^Oin .aY^ C from 
UlJ, the 26Al(n,p)2fMg from 0, the ^^A\{n,af^Nii from O and the 
33S(n,a)30Si from [loj The ^5Cl(n,p)35S rate has been derived from 0, 
the 36Cl(n,p)36S and the ^^C\{n,a)^^F from [Tiijjie 37Ar(n,p)37Cl and the 
3^Ar(n,a)34S from 0. The 39Ar(n,a)36S is from while ^^Ca(n,p)^iK and 
^^Ca(n,a)^^Ar are from [l2^. 



7.4 P decay reactions 



Weak interaction rates (electron captures, f3 and positron decays) are inter- 
polated as a function of the temperature (from 10'' K to 10^° K) and electron 
density (from 1 to 30 cm~^). At temperatures lower than 10^ K we assume a 
constant value e qual to the terrestrial one. For isotopes up to ^^Ar, data have 
been taken from [o^, with the exception of ^Be jl^ and the isomeric state of 
^^Al, for which we refer to i^]- Concerning the unstable isotopes between ^^Ar 
and ^^Ca we use prescriptions by 4^, while between ^^Ca and ^"^Cu (excluding 
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N 



m 
to 



i) w e follow j85|. For ^^Ni and heavier isotopes we use the rates tabulated 
119|, with the exceptions of ''^Se and ^^^Lu, for which we refer respectively 
ijandia. For the few rates not mcluded in any compilation, we use the 



terrestrial value. 



7.5 Isomeric states 



The existence of isomeric s tates of un thermalized isotopes leads to ramifica- 
tions of the s-process flux f jl3lj . jl30l |). In particular, branching points orig- 
inated by the isomeric state of ^^Al, ^^Kr, ^^^Lu and ^^'^Ta require particular 
attention. 

Concerning ^^Al, the proton capture on ^^Mg is split in two distinct reactions: 
the first produces the ground state ^^AF (with terrestrial half life T1/2 = 7.16 
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X 10^ yr) and the second creates the isomer ^e^jm that almost istantaneously 
decays into ^^Mg. 

Concerning ®^Kr, the neutron capture of ^^Kr to ^^Kr™ has a 50% probabihty 
with respect the total cross section (at 30 keV, [IQ]). ^^Kr*" has a non-zero 
probability to decay by internal conversion to its ground state (of 20%), thus 
leading to the following isomeric ratio (IR): 

a(8^Kr(n,7)85Kr™) 

IR = ^ ; ,,,,, ; = 0.42 . 1 

atot(84Kr(n,7)85Kr) ^ ' 

Note that the half life of ^^Kr"^ is 4.48 h, while that of ^^Kr^ is 10.76 yr. 

Concerning ^"^^Lu, below 20 keV we adopt the terrestrial isomeric ratio, namely: 

a(^^^Lu(n,7)^^^Lu-) _ 
^^.o*(^^^Lu(n,7)i^6Lu) 



while above this energy, the partial thermalization through a mediating state 
should also be considered (see 0]) and we assumed an effective IR=0.5. 

Finally, for the complex treatment of the branching between isomeric and 
ground state of ^^°Ta, we refer to jo^ . 



8 Mixing during the third dredge up and formation of the C 
pocket 



In the model presented here we adopt a simple approach based on an argu- 
ment early discussed by P (see also 2^|, 4^|). When the TDU takes place, the 
opacity of the envelope (H-rich) is significantly larger than the opacity of the 
underlying H-exhausted region (He-rich). This fact causes an abrupt change 
of the temperature gradient at the inner border of the penetrating convective 
envelope. In this condition, the convective boundary becomes unstable, be- 
cause any perturbation causing an excess of mixing immediately leads to an 
increase of the opacity and, in turn, to an increase of the temperature gradi- 
ent. This occurrence favours a deeper penetration of the convective instability 
or, in other words, a deeper dredge up. A similar mechanism is responsible 
for the growth of the convective core during central He burning ([o^, jl^l)- A 
different view of the same phenomenon concerns the evaluation of the average 
convective velocity. In the framework of the mixing length theory, this velocity 
is proportional to the difference between the radiative temperature gradient 
(i.e. the gradient necessary to carry out the total energy flux if convection 
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would be inhibited) and the adiabatic temperature gradient. For this reason, 
the average convective velocity usually drops to at the stable boundary of a 
convective layer, where the temperature gradient coincides with the adiabatic 
one. However, when convection penetrates in a region of lower opacity (as 
happens during a third dredge up episode), the difference between the actual 
temperature gradient and the adiabatic gradient grows above and a positive 
average convective velocity is found at the inner border of the convective en- 
velope, which, for this reason, becomes unstable. In principle, as soon as He 
is mixed with the envelope, the opacity and, in turn, the difference between 
the radiative and the adiabatic temperature gradients are reduced. However, 
the mass of the convective envelope is usually much larger than the amount of 
the dredged up material and the relaxing effect of the additional mixing is, in 
practice, negligible (ji^, jl^)- Then, a simple thermodynamic criterion cannot 
be used to determine the real extension of the convective instability. For sure, 
the steep pressure gradient that develops immediately below the formal bor- 
der of the convective envelope limits the penetration of the instability, so that 
the average convective velocity should rapidly drop to 0. In order to mimic 
this behaviour, we assume that in the region underlying the formal convective 
boundary, the average velocity follows an exponential decline, namely: 

(3) 




where d is the distance from the formal convective boundary, Vbce is the velocity 
of the most internal convective mesh. Hp is the pressure scale height and /? is a 
free parameter. Note that this formula is similar to the "overshoot" proposed 
by j61|- However, in our case, since Vbce is usually 0, it produces a negligible 
amount of extra mixing, except during a dredge up. 

Then, the degree of mixing is calculated by means of the following relation: 
=X^ + ^J2,iXk - X°)f,,,AM, (4) 



where the summation is extended over the whole convective zone and the 
superscript ° refers to unmixed abundances. AM^ is the mass of the mesh- 
point k and Mconv is the total mass of the convective zone. The damping factor 
/is: 

kk = ^ (5) 



if At < Tj^k, or 

kk = 1 (6) 
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if At > Tj^k- Here At is the time step and rj^k is the mixing turnover time 
between the mesh-points j and k, namely: 



The mixing velocity (vi) is computed according to the mixing length theory, 
except in the region where the exponential decline is assumed. This algorithm 
allows us to account for the partial mixing that occurs when the time step is 
reduced to or below the mixing timescale. In practice, when the third dredge 
up takes place, complete mixing is obtained within the fully convective zone, 
while the region immediately below, where the convective turnover time scale 
is larger, is only partially mixed. Figure 7 shows the chemical profile, the 
radiative and adiabatic temperature gradients, the average velocity and the 
pressure, in the region around the convective boundary, in a model of 2 Mq, 
during the 5*'* dredge up episode (see section 9). The three major effects of 
the introduction of this exponential decay of the convective velocity are: i) the 
convective boundary is more stable against perturbation; ii) a smooth profile 
of protons within the intershell is left by the TDU; iii) a more efficient TDU 
results. As a consequence of ii), a ^^C pocket almost ^^N-free forms in this 
zone. After a few tests, we found that in order to get a suitable amount of ^^C 
in the pocket, the /3 parameter should be of the order of 0.1. 



9 A thermally pulsing AGB model 

The network and the mixing algorithm described in the previous sections have 
been used to calculate the evolutionary sequence of a 2 Mq star of solar metal- 
licity shown in Figure 2. Before the first thermal pulse, a Reimers' formula for 
the mass loss [r] = 0.4) has been assumed. At this epoch, the total and the 
core masses are 1.95 and 0.55 Mq, respectively. During the thermally pulsing 
AGB phase, we use the mass loss rate prescription discussed in section 2.3 
(see the solid line in Figure 6). The calculation has been stopped when the 
envelope mass was reduced down to ~ 0.3 Mq, the corresponding total mass 
being 0.92 Mq. At this epoch, the TDU ceases (see Figure 8) and hence the 
envelope composition no longer changes. Figure 8 illustrates the evolution of 
the mass coordinates of the maximum energy production within the He in- 
tershell, the maximum energy production within the H burning shell and the 
inner border of the convective envelope. The first TDU episode occurs after 
six TPs, when the core mass is 0.56 Mq. The evolution of 6Mtdv (the amount 
of mass dredged up) and T^^^ (the maximum temperature at the base of the 
convective zone generated by the thermal pulse) are reported in Figure 9. The 
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Fig. 7. Chemical profile {upper panel) , temperature gradients {middle panel), average 
velocity and pressure {lower panel), in the region at the convective boundary during 
the fifth TDU episode. 

possible activation of the ^^Ne(a,n)^^Mg reaction critically depends on T^^x 
that, for this reason, represents a key quantity for the comprehension of the nu- 
cleosynthesis occurring in AGB stars. It is obviously related to the maximum 
He burning luminosity. Then, as for the dredge up, the maximum temper atur e 
is a function of the core mass, the envelope mass and the composition lisl . 



As expected, the deepness of the third dredge up initially increases, due to the 
increase of the core mass. Then, during the TP-AGB phase, the effect of the 
envelope erosion becomes important, and the dredge up efficiency decreases, 
dropping to after the last computed thermal pulse. C/0 becomes larger than 
1 when the core mass is about 0.6 Mq and the luminosity is log L/Lq = 3.85, 
which corresponds to a bolometric magnitude of —4.9 mag. We recall that the 
majority of the C(N) stars in the galactic disk have bolometric magnitudes 
around —5 mag. The final C/0 is 1.67. 

'^max initially increases and attains a nearly asymptotic upper limit, of about 
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Fig. 8. The evolution throughout the TP-AGB phase of a star with initial mass M 
= 2 Mq and solar metallicity. The three lines represent (from the top): the inner 
border of the convectivc envelope, the location of the maximum energy production 
within the H burning shell and the location of the maximum energy production 
within the He-rich intershell. 

285 X 10^ K. At such a low temperature the ^^Ne(Q;,n)^^Mg is marginally ac- 
tivated. On the contrary, an important s-process nucleosynthesis is driven by 
the ^^C(a,n)^^0 reaction. A ^^C pocket forms after each TDU episode. The 
^^C in the first two pockets is only partially burned during the interpulse and 
the residual is engulfed into the convective zone generated by the subsequent 
thermal pulse. In all other cases, the ^^C is fully consumed during the inter- 
pulse. The outward shift of the location of the maximum energy production 
within the intershell, in Figure 8, marks the onset of the ^^C burning during 
the interpulse. Figure 10 illustrates the various steps of the formation of the 
^th i3q pocket. Note that the zone where the ^^C pocket forms (panel d)) is 
partially overlapped by a more external thin ^^N pocket. The maximum neu- 
tron density (about 10^ cm~^) is attained in the more internal layer of the ^^C 
pocket, where the -"^^N is less abundant. The effective abundances of ^^C in 
all the pockets we found (defined as X^s{^^C)=X{^^C) - X{^^N) x 13/14) are 
reported in Figure 11. The extension of the pockets decreases with time, the 
first one beeing the largest. 

The ashes of the neutron-capture nucleosynthesis allowed by the ^^C burning 
are spread within the intershell by the subsequent convective thermal pulse. 
Later on, as a consequence of the third dredge up, the envelope composition 
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Fig. 9. Evolution of the maximum temperature at the base of the convective zone 
generated by the thermal pulse (panel a) and dredged up mass per pulse (panel b) 
during the AGB phase of the model of Figure 1. 

is polluted with the products of the s process. The final surface composition 
resulting after nine TDU episodes is reported in Figure 12. Among the light 
elements, note the logarithmic enhancements of C (0.5 dex), F (0.3 dex), 
Ne (0.14 dex) and Na (0.12 dex). We recall that C is the main product of 
the 3Q;-burning. Fluorine is synthesized in the convective shell generated by 
the thermal pulse via the ^^N(q;,7)^^F reaction, the ^^N being accumulated 
in the -"^'^C pocket mainly via the chain -'^^N(n,p)^^C(«,j9)"^^0(p, a)"^^N^ . The 
nitrogen enhancement in the envelope (0.3 dex) is totally due to the first 
dredge up. In fact, the ^^N left by the H burning within the intershell is fully 
converted into ^^Ne during the convective thermal pulse, giving rise to the 
resulting Ne enhancement. The Na enhancement is partly due to the proton 
captures on ^^Ne, occurring within the H-burning shell, and partly due to the 



' A minor contribution to the synthesis of the in the pocket may eventually 
come from the ^^N(n,7)^^N and from the ^^N(j», 7)^^0 reactions. 



23 




Fig. 10. The variation of the relevant chemical species in the region where the 
^^C pocket forms after the occurrence of the third dredge up. Each panel refers to 
different time steps, namely: At = 78, 1000, 2000 and 14000 years, respectively. 



neutron captures on ^^Ne followed by the ^^Ne decay, taking place within the 
He intershell. Concerning the s-process nucleosynthesis, all the elements (from 
Sr to Pb) are enhanced at various degree. We find that the abundance of Sr, 
Y and Zr at the first s peak, the so-called Is elements {light s elements), is 
comparable with the one of the hs elements [heavy s elements) Ba, La, Ce, 
Pr, Nd at the second s peak. Lead is underproduced with respect to barium, 
as expected for AGB stars of this metallicity (see next section). In Figure 
13, the predicted [hs/ls] when (C/0=1) is compared to those measured in 
galactic C(N) giants Note that, for a given metallicity, this intrinsic index 
is indicative of the physical conditions of the s-process site. 
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Fig. 11. The effective (see text) mass fraction of ^^C within all the pockets is re- 
ported. Each pocket has been shifted in mass in order to superimpose their external 
borders. The point of the abscissa is arbitrary. The extension of the pocket de- 
creases with time. The first is about 6.5 • 10~^ Mq large. The last is about 6 time 
smaller. 




Z 



Fig. 12. Final surface composition of the model in Figure 1. The standard spectro- 
scopic notation, namely [El/Fe]=log(A^(El)/iV(Fe))-log(iV(El)/Af(Fe))o, is used, 
where El is the name of a generic element. 
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Fig. 13. Observed [hs/ls] ratio of a sample of galactic C(N) stars, compared with our 
theoretical prediction. The huge cross represents typical observational error bars of 
C(N) stars, whose spectra are very crowded (see 

10 Nucleosynthesis in low-metallicity AGB stars: lead stars 



The last challenging chapter of the AGB nucleosynthesis story is provided by 
the high resolution spectroscopy of very metal-poor C- and s-rich stars. The 
present generation of halo stars is old (~ 14 Gyr) and, therefore, is made of low 
mass objects (M < 0.9 Mq). Then, when a low mass star reaches the AGB, 
its envelope is so small that the TDU never takes place. However, available 
spectroscopic surveys for very metal-poor stars 0) found that ~20% 

to 30% of the candidates ([Fe/H] > —2.5) are carbon rich. These stars are 
probably low mass dwarfs or giants, with a lifetime comparable with the age 
of the Galaxy, belonging to binary systems. It is possible that in most cases the 
C enhancement is the result of an ancient accretion process (by stellar wind or 
Roche lobe overflow j?^) from a more massive AGB companion (now a white 
dwarf). In such a case, also the products of the neutron capture nucleosynthesis 
were enhanced in the material accreted onto the secondary star. 

In low mass AGBs, the main source of neutrons, the ^'^C(Q;,n)^^0, is primary- 
like, i.e. not directly affected by the metallicity of the pristine material. Nev- 
ertheless, the iron seeds scale with the metallicity, so that the lower the metal- 
licity the larger is the number of neutrons available per seed. As a result, at 
low metallicities, most of the seeds are converted into ^°^Pb, at the termina- 
tion point of the s-process fiuence. Note that in stars suffering recurrent TDU 
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Fig. 14. Predicted [hs/ls] versus [Fe/H] for different ^^C pocket choices. 

episodes, the amount of ^^C in the envelope increases with time. This primary 
carbon is mainly converted into ^^N, within the H-burning shell, and, later 
on, into ^^Ne, within the convective zone generated by a thermal pulse (see 
section 6). In very metal poor stars, the primary neon accumulated in the in- 
tershell becomes the main seed for the s process (see also IHsl). For these 
reasons, in a very metal-poor AGB star experiencing a few TDU episodes, a 
consistent enhancement of lead is expected High-resolution spectroscopy 
of very metal-poor C-rich stars largely confirms such a qualitative expecta- 
tion. However, the precise prediction of the heavy element enhancements at 
the surface of the secondary star is more complex than previously sketched. 
We have to recall that we are looking at the intershell material that was mixed 
with the envelope of the primary AGB stars during the various TDU episodes 
and, later on, further diluted within the pristine material of the envelope of 
the secondary star. If this secondary star (actually the C-rich object presently 
detected) , is evolved off the main sequence, this late dilution may be particu- 
larly efficient. In addition, the amount of mass accreted depends on the orbital 
parameters of the system that are, in most cases, unknown. 



In this context, in order to investigate the efficiency of the s-p rocess site, 
intrinsic indices [hs/ls] and [Pb/hs] are particularly useful ( j50| , [s^). The 
variation of these two indices with metallicity are shown in Figures 14 and 15. 
These calculations are based on a grid of older FRANEC models, where the 
Reimers's parameterisation [rj = 0.5) was used and no velocity profile below 
the convective envelope was introduced. The n-capture nucleosynthesis was 
calculated with a post-process code that follows in detail the physical and 
thermodynamical structure of the He-rich intershell in the TP-AGB phase. 
After each third dredge up, diffusion of protons in a thin layer in the top re- 
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gion of the He-rich intershell was assumed. It gives rise to a ^^C pocket at H 
reignition; the amount of ^^C in the pocket and its profile in mass was taken 
as a free parameter. The standard case (ST, see ji^) was shown to reproduce 
the main s-process component in the solar system for low mass AGB stars and 
half-solar metallicity. The ^^C pocket efficiency was assumed identical for all 
thermal pulses followed by TDU. Neutrons released by the ^^C pocket are of 
primary origin, depending on proton captures on the freshly nucleosynthesised 
^^C. Varying the metallicity, the s-process distribution changes because it de- 
pends on the number of neutron captured per seed and, at low metallicities, 
peaks at Pb, at the termination of the s-process path. As recalled in section 
7.2, a large spread of ^^C pocket efficiencies is needed in order to explain 
the spectroscopic observations of MS-S-C(N) stars belonging to the galactic 
disk In Figures 14, 15, different lines correspond to different ^^C pockets, 
where the ST case was multiplied or divided by a given factor as indicated 
in the inset. The STx2 case is an upper limit because ingestion of a larger 
abundance of protons would favour the ^^N production by proton capture on 
^^C. As a comparison, spectroscopic measurements of C and s-rich metal-poor 
stars are reported in Figure 15. It results that for very metal-poor stars, [hs/ls] 
is weakly affected by a change of the efficiency of the ^^C pocket, while large 
spreads of [Pb/hs] are expected, in quite satisfactory agreement with observa- 
tions. This is in accord with the large spread of ^^C pocket efficiencies required 
by the observation of MS-S-C(N) stars (Q, Q and 0). 



Q. 




ST*2 
ST*1 .3 
ST 
ST/1 .5 
ST/2 
ST/3 
ST/6 
ST/12 
ST/24 
ST/45 
ST/75 



M= 1.5M„ 



-3 -2.5 -2 -1.5 -1 -0.5 0.5 

[Fe/H] 

Fig. 15. Predicted [Pb/hs] versus [Fe/H] for different ^^C pocket choices. Spectro- 
scopic data of s-rich and lead-rich stars are included for comparison (adapted from 



Very low metallicity s-rich stars have an additional peculiarity: they are usually 
also N-rich. The enhancement of nitrogen could be due to the hot bottom 
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burning, but this would exclude low mass AGB stars that have a too cool 
temperature at the base of the convective envelope. The activation of the 
so-called cool bottom process during the AGB [91j, a slow deep circulation 
taking place in the radiative region located between the H-burning shell and 
the inner border of the convective envelope, could contribute to the systematic 
occurrence of the nitrogen overabundance. More likely, [g^ (see also |47| ) 
suggested that the scarcity of CNO in very metal-poor stars may allow the 
convective zone generated by the first strong TP to ingest some protons from 
the envelope. In this hot CN cycle would take place within the He-rich 

intershell and consistent amounts of ^^C and ^^N are produced. More recently 
[tH - on the basis of detailed stellar models computations, proposed that a 
peculiar s-process nucleosynthesis, characterized by a high neutron density, 
could be activated by a-capture on the freshly synthesized ^^C. The detailed 
s-process nucleosynthes is fo r a low-mass and very metal-poor AGB model 
has been calculated by 114| . We found that the ingestion of proton into the 
convective intershell region and the subsequent deep TDU provide the required 
mechanism to explain the N enhancement, but the high neutron density is 
maintained only for a very short time (few days), with negligible consequences 
on the overall surface overabundances. 



11 Conclusion. 



We have reviewed the theoretical investigation of nucleosynthesis and evolu- 
tion of low mass AGB stars. By introducing a full network (from H to Bi) in 
our stellar evolution code, we have obtained a self-consistent nucleosynthesis 
calculation for a 2 Mq stellar model with solar composition. 

Below the convective envelope of a thermally pulsing AGB star, the H burning, 
the He burning and the s process modify the internal composition. Then, the 
occurrence of recurrent third dredge up episodes induces the modification of 
the surface composition. The adoption of an exponential decay of the average 
convective velocity at the inner border of the convective envelope allows us to 
follow the formation of a ^'^C pocket in the top layer of the He-rich intershell. 
During the interpulse, this pocket is compressed and heated, untill neutrons 
are released via the ^^C(a,n)^^0 reactions. This is the main neutron source, 
responsible for the s-process nucleosynthesis in low mass AGB stars. A second 
(marginal) neutron burst, via the ^^Ne(a,n)^^Mg reactions, may eventually 
occurs at the base of the convective zone generated by a thermal pulse. One of 
the main source of uncertainty in modeling AGB stars concerns the choice of 
the mass loss rate. A new semiempirical calibration of the mass loss has been 
derived from a collection of the most recent observational data of variable AGB 
stars. This calibration has been used for the AGB model calculation presented 
here. We have followed the evolution up to the last third dredge up episode. 
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The final elemental distribution is representative of the one expected for the 
intrinsic carbon stars observed in the disk of the Milky Way. A comparison 
with available spectroscopic analysis shows a reasonable agreement. 

40 years h ave p assed since the discovery of the thermal pulses by Schwarzschild 
& Harm |l09j and many efforts have been made for the comprehension of 
the AGB evolution and nucleosynthesis. The progressive improvement of the 
computational capability allows the development of more complex theoretical 
tools. The main tasks for the next future will be the comprehension of the role 
played by hydrodynamical phenomena usu ally neglected in standard stellar 
model computations (see e.g H, H, liij, 0, 0, 0) and the refinement 
of the nuclear physics input data. 

We wish to express our thanks to Debora Delaude and Marco Pignatari for 
helping in the preparation of this manuscript. We are indebted to Franz 
Kappeler and Inma Dominguez for a long profitable and continuous interac- 
tion. Part of this work has been supported by the Italian MIUR-FIRB Project 
"The Astrophysical Origin of the Heavy Elements beyond Iron" and by the 
Italian MIUR-PRIN Project 2004 "Nuclear astrophysics in low mass stars". 
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